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The Observed Growth of Massive Galaxy Clusters IV: 
Robust Constraints on Neutrino Properties 



o 
a 

(N 

o 

u 

Or 
6 

o3 



> 
oo 
oo 



ON 

o 



13 



A. Mantz, 1 ' 2 ' 3 * S. W. Allen 2 ' 3 and D. Rapetti 2 ' 3 

1 NASA Goddard Space Flight Center, Code 662, Greenbelt, MD 20771, USA 

2 Kavli Institute for Particle Astrophysics and Cosmology, Stanford University, 452 Lomita Mall, Stanford, CA 94305-4085, USA 
3 SLAC National Accelerator Laboratory, 2575 Sand Hill Road, Menlo Park, CA 94025, USA 



Accepted 2010 April 6. Received 2010 March 14; in original form 2009 November 5 



ABSTRACT 

This is the fourth of a series of papers in which we derive simultaneous constraints on 
cosmological parameters and X-ray scaling relations using observations of the growth of 
massive, X-ray flux-selected galaxy clusters. Our data set consists of 238 clusters drawn 
from the R OSA T All- Sky Survey, and incorporates extensive follow-up observations 
using the Chandra X-ray Observatory Here we examine the constraints on neutrino 
properties that are enabled by the precise and robust constraint on the amplitude of 
the matter power spectrum at low redshift available from our data. In combination 
with cluster gas mass fraction, cosmic microwave background, supernova and baryon 
acoustic oscillation data, and incorporating conservative allowances for systematic 
uncertainties, we limit the species-summed neutrino mass, M v , to < 0.33 eV at 95.4 
per cent confidence in a spatially flat, cosmological constant (ACDM) model. In a flat 
ACDM model where the effective number of neutrino species, N e g, is allowed to vary, 
we find N e g = 3.4+q'j (68.3 per cent confidence, incorporating a direct constraint on 
the Hubble parameter from Cepheid and supernova data). We also obtain results with 
additional degrees of freedom in the cosmological model, in the form of global spatial 
curvature (Ok) an d a primordial spectrum of tensor perturbations (r and n t ). The 
results arc not immune to these generalizations; however, in the most general case we 
consider, in which M„, N e g, curvature and tensors are all free, we still obtain M v < 
0.70 eV and N g = 3.7 ± 0.7 (at, respectively, the same confidence levels as above). 
These results agree well with recent work using independent data, and highlight the 
importance of measuring cosmic structure and expansion at low as well as high (z ~ 
1100) redshifts. Although our cluster data extend to redshift z = 0.5, the direct effect 
of neutrino mass on the growth of structure at late times is not yet detected at a 
significant level. 
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1 INTRODUCTION 

Observations of neutrino flavor oscillation have con- 
clusively shown that the neutrino mass eigenstates 
are non-degenerate (e.g. Fukuda et al. 1998, 2002; 
Ahmad et al. 2002; Ahn et al. 2003, 2006; Eguchi et al. 
2003; Sanchez et al. 2003; Giacomelli & Margiotta 2004; 
Aharmim et al. 2005). Although, these observations can 
place tight constraints on the differences in squared 
mass, measuring the absolute mass scale remains chal- 
lenging. Current laboratory efforts focus on tritium beta 
decay (e.g. Lobashev 2003; Kraus et al. 2005) and neu- 



trinoless double beta decay (e.g. Aalseth et al. 1999; 
Klapdor-Kleingrothaus et al. 2001; Arnaboldi et al. 2005; 
Arnold et al. 2005); the latter, if observed, would addi- 
tionally indicate that neutrinos are Majorana rather than 
Dirac fermions. The squared mass differences measured 
from flavor oscillations place a lower bound on the sum of 
the three masses, M u = m,j, at ~ 0.056 (0.095) eV/c 2 in 
the normal (inverted) hierarchy, while current tritium beta 
decay results provide an upper bound on the mass of the 
electron neutrino at ~ 2eV/c 2 (thus on M v at ~ 6 6V/C 2 ). 1 
Because neutrinos play a prominent role in the early 
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1 Henceforth, we set c = 1. 
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Universe, cosmological observations are also sensitive to 
their properties (for a review, see Lesgourgues & Pastor 
2006; see also Section 2). The primary effect of non-zero 
neutrino mass on cosmological observables is to suppress 
the formation of cosmic structure on intermediate and small 
scales. Comparison of the Cosmic Microwave Background 
(CMB), which reflects large-scale structure at early times, 
with measurements of the intermediate- or small-scale struc- 
ture in the local Universe thus provides a way to constrain 
the absolute mass scale of the neutrino (e.g. Fukugita et al. 
2000). 

This approach has the disadvantage that any neutrino 
properties inferred are at some level sensitive to our incom- 
plete understanding of cosmology, in particular dark energy 
and inflation. It has long been recognized, however, that us- 
ing complementary measurements of structure, as described 
above, significantly reduces the sensitivity of the results to 
such necessary assumptions (e.g. Allen, Schmidt, & Bridle 
2003; Tegmark et al. 2004). Nevertheless, it is imperative 
that any systematic uncertainties affecting the measure- 
ments of cosmic structure be properly accounted for in order 
to obtain robust results. 

The amount of structure in the local Universe is 
generally described by the parameter as, defined as the 
present day root-mean-square fluctuation of the linearly 
evolved matter density field, smoothed by a spherical 
top-hat window of comoving radius 8/i _1 Mpc; here h = 
-ffo/100 km s" 1 Mpc -1 is the normalized Hubble parame- 
ter. At present, the most robust observation for measur- 
ing this quantity is arguably the abundance of massive 
galaxy clusters; 2 recent advances in cluster simulation (e.g. 
Nagai, Vikhlinin, & Kravtsov 2007), comparisons of differ- 
ent mass measurement techniques (e.g. Bradac et al. 2008; 
Mahdavi et al. 2008; Newman et al. 2009), and the avail- 
ability of robust mass proxies (e.g. Allen et al. 2004, 2008; 
Kravtsov, Vikhlinin, & Nagai 2006) have significantly re- 
duced systematic uncertainties associated with the mea- 
surement of cluster masses. Consequently, recent estimates 
of as based on independent analyses of galaxy cluster 
data (both optically and X-ray selected clusters) are in 
very good agreement [a s ~ 0.8; Mantz et al. 2008, 2009b 
(hereafter Paper I); Henry et al. 2009; Vikhlinin et al. 2009; 
Rozo et al. 2010]. Recent analyses of cosmic shear data 
are also in good agreement (Benjamin et al. 2007; Fu et al. 
2008), and provide compatible constraints on the neu- 
trino mass to our own (Tereno et al. 2009). Galaxy 
redshift surveys have also produced comparable results 
(Thomas, Abdalla, & Lahav 2009). 

The number of neutrino species participating in weak 
interactions is known to be three to high precision 
(Amsler et al. 2008). However, the possibility remains that 
additional "sterile" species exist. Results consistent with the 
existence of a fourth neutrino were reported by the Liquid 
Scintillator Neutrino Detector (LSND; Aguilar et al. 2001); 
these results are disfavored by the MiniBooNE experiment 



2 We note that current cluster measurements do not constrain erg 
independent of spectral index of the power spectrum, n B . How- 
ever, the best fitting value of as does not vary rapidly with n B ; 
moreover, n s is well constrained by CMB data in all cosmological 
models considered here. 



( Aguilar- Arevalo et al. 2009a, b), although the interpreta- 
tion remains somewhat ambiguous at present. This is an- 
other question that cosmological observations can address. 
In particular, the synthesis of light elements is sensitive to 
the number of relativistic species present in the early Uni- 
verse, since these determine the expansion rate at that time; 
however, observations of primordial deuterium and helium 
abundances are challenging and are subject to large system- 
atic uncertainties. An independent probe is provided by the 
CMB, in combination with other cosmological data, as de- 
scribed in Section 2. In this work, we consider only the case 
where the neutrino species (whatever their number) have 
approximately degenerate mass; in this case, M v and the 
effective number of neutrino species, N c g, are sufficient to 
describe the cosmological effect of neutrinos. More general 
mass splittings, and in particular the case favored by the 
initial LSND results, in which a sterile neutrino is signifi- 
cantly more massive than the other species, require a more 
complete treatment (as in, e.g., Crotty et al. 2004). 

In this paper, we apply the statistically rigorous analysis 
method of Paper I and the X-ray flux-limited cluster sam- 
ples and follow-up observations described in Mantz et al. 
(2009a, hereafter Paper II) to the problem of inferring neu- 
trino properties from cosmological data. These data (collec- 
tively referred to here as the cluster X-ray Luminosity Func- 
tion, or XLF) provide a robust means to measure as', our 
analysis method includes generous systematic allowances 
and accounts fully for all parameter degeneracies. In ob- 
taining our results, we also incorporate CMB data and mea- 
surements of cosmic distance in the form of cluster gas mass 
fractions (/ ga s), type la supernova (SNIa) fluxes and Baryon 
Acoustic Oscillation (BAO) data (see Section 3 and refer- 
ences therein). We note that Reid et al. (2010) obtain very 
similar results to our own by importance sampling 5-year 
Wilkinson Microwave Anisotropy Probe ( WMAP) results us- 
ing a prior based on the analysis of optically selected clusters 
by Rozo et al. (2010). 

The basic cosmology that we consider in this paper 
is the spatially flat, cosmological constant (ACDM) model 
parametrized by the mean baryon density, Qb', the mean 
total matter density including baryons, neutrinos and cold 
dark matter (CDM), f2 m ; the Hubble parameter, Ho; the 
normalization of the matter power spectrum, as; the spec- 
tral index of the primordial scalar power spectrum, n s ; and 
the optical depth to reionization, r. Here the mean densities 
refer to the present day (redshift z = 0), since their values at 
other times are then determined by the Friedmann equation, 
and as is the 2 = variance in the matter density field at 
scales of 8/i _1 Mpc, as defined above. This simple and com- 
monly used model assumes M v = and N c s = 3.046, the 
predicted value for the three weakly interacting neutrinos. 

In addition to freeing M„ and Nes, we will consider 
the effect of marginalizing over other parameters with which 
they are degenerate. However, since the flat ACDM model 
is known to provide a good fit to currently available cos- 
mological data, we are conservative in incorporating these 
additional degrees of freedom. In particular, we consider gen- 
eralizing the description of dark energy, either by allowing 
global spatial curvature or by varying the dark energy equa- 
tion of state in a flat universe, and including the effects of 
primordial tensor perturbations. The former case is param- 
eterized by the effective curvature energy density, fik , or the 
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equation of state parameter, w, while in the latter case we 
simultaneously marginalize over the tensor-to-scalar ratio, r 
(defined with respect to wavenumber k = 0.002ft Mpc -1 , as 
in, e.g., Komatsu et al. 2009), and the tensor power spectral 
index, n t . 

In our results, we will consistently quote one-sided lim- 
its (upper bounds) on M v and r at the 95.4 per cent confi- 
dence level and two-sided constraints on all other parameters 
at 68.3 per cent confidence. 



2 BACKGROUND 

Although the purpose of this work is to explore the con- 
straints enabled by the XLF data, it is instructive to begin 
by considering the cosmological constraints on M v that have 
been obtained without using measurements of local cosmic 
structure. As described by Ichikawa, Fukugita, & Kawasaki 
(2005), Dunkley et al. (2009) and Komatsu et al. (2009), 
CMB observations alone provide a limit that, while weak 
(M„ < 1.3 eV at 95.4 per cent confidence), is relatively ro- 
bust to assumptions about dark energy and primordial ten- 
sors. 3 While we do not review the details here, the constraint 
arises because the effect of neutrino mass on the CMB tem- 
perature anisotropy spectrum is qualitatively different when 
neutrinos are relativistic compared with non-relativistic dur- 
ing the decoupling epoch. In particular, for M v < 1.5 eV, 
changes in M v can be easily mimicked by a corresponding 
change in the Hubble parameter, Ho- As a result, M u is 
strongly degenerate with Ho in analyses using only CMB 
data, which does not itself constrain Ho- The combination of 
CMB data with distance measurements in the form of SNIa 
and BAO data can place a constraint on Ho, improving the 
limits on M„ by roughly a factor of two (M„ < 0.67 eV for 
a spatially flat ACDM background; Komatsu et al. 2009). 

As we will show, the inclusion of a measurement of lo- 
cal cosmic structure in the form of galaxy cluster abundance 
can improve these limits by another factor of two or greater, 
and furthermore greatly increase the robustness of the re- 
sults to assumptions about the nature of dark energy and 
inflation. This improvement is possible because the combina- 
tion of large-scale, high-redshift and intermediate-scale, low- 
redshift measurements of cosmic structure exploits the defin- 
ing characteristic of light, weakly interacting particles: they 
are relativistic during the earliest stages of cosmic structure 
formation but non-relativistic at the present day. 4 The im- 
print of neutrinos on cosmic structure, and in particular on 
the growth of structure from the surface of last scattering 
to the present, thus differs from that of the photon back- 
ground (relativistic at all times) and CDM (assumed to be 
non-relativistic at all times). In particular, a non-zero mass 
in neutrinos results in a net suppression of the growth of 
structure on scales smaller than their free-streaming length 

3 We note that this constraint is not entirely robust, and can be 
significantly degraded when multiple additional degrees of free- 
dom are included in the cosmological model. In several of the 
models we consider, the combination of CMB data and cosmic 
distance measurements produces weaker constraints; however, in 
the worst case, the upper limit is still ~ 2.5 eV (Tables 1 and 2). 

4 The latter condition follows from the lower bound M„ > 
0.056 eV obtained from flavor oscillation measurements. 



(approximately corresponding to wavenumber 10 _2 ftMpc _1 
today) relative to an equal mass in CDM (Bond et al. 1980). 

In practical terms, the constraint on M„ arising from 
the combination of CMB and cluster data can be under- 
stood as follows. Even when M v =^ 0, the CMB data pro- 
vide a good constraint on the high-redshift amplitude of the 
power spectrum, since the direct effect of neutrino mass on 
the observed temperature anisotropies is relatively small. 
However, the late-time, intermediate-scale power spectrum 
amplitude (the value of as) that would be predicted based on 
that high-redshift constraint is very sensitive to M v , since 
neutrino mass suppresses the development of structure on 
those scales. This strong degeneracy is apparent in the blue 
confidence regions in Figure 1. By constraining as indepen- 
dently of M v , clusters (or other low-redshift structure mea- 
surements) break this degeneracy, improving the constraints 
on neutrino mass (gold contours in the figure). 

Even more power is available if the growth of structure 
can be measured, exploiting the time-dependent nature of 
the effects of neutrino mass on cosmic structure develop- 
ment (e.g. Lesgourgues & Pastor 2006). However, we show 
in Section 5.1 that the growth of structure in current cluster 
data, while able to constrain models of dark energy (Paper I; 
Vikhlinin et al. 2009) and modified gravity (Rapetti et al. 
2009a,b; Schmidt, Vikhlinin, & Hu 2009), are not yet suffi- 
cient to contribute to constraints on M v . 

In the case of the effective number of relativistic species, 
N e a, CMB data alone can provide a lower bound, due to the 
gravitational effect of their anisotropic stress. Apart from 
providing this bound, however, the CMB data are ham- 
pered by a near-perfect degeneracy between N e ff and Qmh 2 , 
the physical matter density (Dunkley et al. 2009). This de- 
generacy arises because CMB observations do not constrain 
either N c g or f2 m ft 2 directly, only the redshift of matter- 
radiation equality, z cq , which is a function of both param- 
eters. Again, the addition of cosmic distance measurements 
can reduce the effect of this degeneracy by placing an in- 
dependent constraint on fi m ft 2 ; Komatsu et al. (2009) ob- 
tained iVcff = 4.4 ±1.5 (68 per cent confidence). 

Measurements of local structure can also contribute 
to improved constraints on A^ff, as we show in Section 6. 
The contribution of clusters to this problem is less direct 
than in the case of M v : the constraint on as and f2 m from 
cluster data improves the constraint on f2 m ft 2 from the 
CMB+distance measurements (e.g. Figure 1 of Paper I), in- 
directly improving the constraint on N e g through the de- 
generacy described above. 



3 DATA 

The galaxy cluster data used in this work, as well 
as their selection and reduction, are discussed in de- 
tail in Paper II. Using three wide-area cluster samples 
drawn from the ROSAT All-Sky Survey (RASS; Triimper 
1993) - the ROSAT Brightest Cluster Sample (BCS; 
Ebeling et al. 1998), the ROSAT-ESO Flux-Limited X-ray 
sample (REFLEX; Bohringer et al. 2004), and the bright 
sub-sample of the Massive Cluster Survey (Bright MACS; 
Ebeling, Edge, & Henry 2001; Ebeling et al. 2010) - we se- 
lect a statistically complete sample of 238 X-ray luminous 
clusters covering the redshift range z < 0.5. Of these 238 
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clusters, 94 have follow-up Chandra or ROSAT observations 
that we incorporate into the analysis. From the follow-up 
observations, we measure X-ray luminosity, average temper- 
ature, and gas mass within rsoo. 5 The gas mass is used as a 
proxy for total mass, using the finding of Allen et al. (2008) 
that the gas mass fraction, / gas = M gas /M t ot, is a constant 
for hot, massive clusters (see also Section 4). Cluster / gas 
measurements additionally provide a precise measure of cos- 
mic distance, so we include the full data set and analysis of 
Allen et al. in this work. 

In addition to these cluster data sets, we incorporate 
CMB, SNIa and BAO data. Our analysis of the CMB 
anisotropics uses 5-year WMAP data (Hinshaw et al. 2009; 
Hill et al. 2009; Nolta et al. 2009) with the March 2008 ver- 
sion of the WMAP likelihood code 6 (Dunkley et al. 2009), 
as well as Arcminute Cosmology Bolometer Array Receiver 
(ACBAR) data at smaller angular scales (Reichardt et al. 
2009) and measurements of the CMB polarization from the 
Background Imaging of Cosmic Extragalactic Polarization 
(BICEP) instrument (Chiang et al. 2010). The SNIa results 
are derived from the Union compilation (Kowalski et al. 
2008), which includes data from a variety of sources (307 
SNIa in total; see references in Paper I). Our analysis of 
BAO data uses the constraints on the ratio of the sound 
horizon to the distance scale at z = 0.25 and z = 0.35 
derived by Percival et al. (2007) from the galaxy correla- 
tion function in 2dF (Colless et al. 2001, 2003) and Sloan 
Digital Sky Survey (Adelman-McCarthy et al. 2007) data. 
Finally, we employ in some cases a Gaussian prior on the 
Hubble parameter, h — 0.742 ± 0.036, based on the results 
of Riess et al. (2009). 7 



4 ANALYSIS METHOD AND SYSTEM ATICS 

Our results are obtained via Markov Chain Monte Carlo 
(MCMC), employing the Metropolis sampler embedded in 
the COSMOMC code of Lewis & Bridle (2002). 8 The May 2008 
COSMOMC release includes the 5-year WMAP and Union su- 
pernova data and analysis codes; an additional module im- 
plementing the /g as analysis has also been publicly released 
(Rapetti, Allen, & Weller 2005; Allen et al. 2008). 9 Further 
modifications were made to include the likelihood codes for 
the XLF and BAO data. The CMB and matter power spec- 
trum calculations were performed using the CAMB package 
of Lewis, Challinor, & Lasenby (2000). 10 

When analyzing the CMB data, we marginalize over a 
plausible range in the amplitude of the Sunyaev-Zel'dovich 
signal due to galaxy clusters (0 < Asz < 2; introduced by 
Spergel et al. 2007). Our analysis of the Union supernova 



rsoo is defined to be the radius within which the mean enclosed 
density is 500 times the critical density at the cluster's redshift. 

6 http://lambda.gsfc.nasa.gov 

7 Strictly speaking, this result should be interpreted as a con- 
straint on the luminosity distance to z ~ 0.04, as noted by 
Reid et al. (2010). However, as those authors conclude, the dis- 
tinction between that approach and a straightforward prior on h 
is very small in practice. 

8 http : //cosmologist . info/ cosmomc/ 

9 http : / /www. Stanford . edu/~ drapetti/f gas_module/ 

10 http://www.camb.info 



sample of Kowalski et al. (2008) includes their treatment of 
systematic uncertainties, which accounts for the effects of 
Malmquist bias and uncertainties in lightcurve fitting and 
photometry (among others). 

The method used to analyze cluster / gas data is de- 
scribed in full by Allen et al. (2008). It incorporates gen- 
erous systematic allowances for instrument calibration (10 
per cent), non-thermal pressure support (10 per cent, 
Nagai et al. 2007), the depletion of baryons in clusters with 
respect to the cosmic mean (20 per cent), and evolution with 
redshift of the baryonic and stellar content of clusters (10 
and 20 per cent). 

The analysis of the XLF data is detailed in Paper I. The 
method combines cluster survey data with follow-up obser- 
vations in an internally consistent way, rigorously account- 
ing for the effects of Malmquist and Eddington biases and 
parameter degeneracies. Conservative systematic allowances 
are included to account for uncertainty in the predicted clus- 
ter mass function, the overall cluster survey completeness 
and purity, and instrument calibration. 

As discussed in Section 2 and illustrated in Sections 5 
and 6, the XLF data contribute to constraints on M v and 
Aefj through their measurement of as- The posterior uncer- 
tainty on this measurement (from the XLF alone) is ~ 6 
per cent, and is dominated by systematic uncertainty, as 
described in Paper I. Specifically, the uncertainty in erg is 
determined by the accuracy with which cluster masses can 
be measured. As detailed in Paper II, we do not directly 
infer cluster masses at rsoo by assuming hydrostatic equilib- 
rium of the intracluster medium, a procedure which, when 
applied to a typical cluster, introduces a large and highly 
variable bias (Faltenbacher et al. 2005; Rasia et al. 2006; 
Nagai et al. 2007). Instead, we use the gas mass at rsoo, 
which can be measured without significant bias, as a proxy 
for total mass. The proxy relation is calibrated using the 
/g as data of Allen et al. (2008), which effectively consist of 
gas mass and total mass measurements for the subset of 
clusters where the hydrostatic assumption is applicable (the 
most massive, dynamically relaxed clusters). Of the system- 
atic allowances described above, the most relevant for this 
procedure are the allowances for non-thermal support and 
instrument calibration in the / gas analysis. In addition, we 
account for uncertainty in the difference in / gas between 
?~25oo ( as measured by Allen et al.) and rsoo, as well as pos- 
sible scatter in / gaB from cluster to cluster (see Paper II). 
Ultimately, both our estimates of individual masses and the 
mean cluster mass scale include a systematic error budget 
of ~ 15 per cent. 



5 LIMITS ON M v 
5.1 Simple models 

We first consider the case of a spatially fiat ACDM cosmol- 
ogy with non-zero neutrino mass. For this model, the joint 
68.3 and 95.4 per cent confidence regions in the M v -o% plane 
from the combination of CMB, / gas , SNIa and BAO data ap- 
pear in the left panel of Figure 1 (blue contours). The 95.4 
per cent confidence upper bound on M v from these data (in- 
cluding the systematic allowances described in Section 4) is 
M v < 0.61 eV (Table 1), a marginal improvement over that 
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Sm„ (eV) Sm„ (eV) 

Figure 1. Joint 68.3 and 95.4 per cent confidence regions in the M v -as plane from the combination of CMB, / ga s, SNIa and BAO data 
(blue), and the combination of those data with the XLF (gold). The XLF data provide a tight constraint on as, breaking the degeneracy 
in this plane. To demonstrate the robustness of constraints on M v obtained using the XLF data, we compare (left panel) results for a 
simple ACDM+M„ model with (right panel) results obtained when nuisance parameters are included in the model (in this case, Q^, r 
and Tit; note the difference in scale). Note that conservative systematic uncertainties are included here (and in all subsequent figures). 



obtained by Komatsu et al. (2009) using only WMAP, SNIa 
and BAO data. 

However, it is clear in the figure that a tight constraint 
on as can improve the limits on M v . The upper limit ob- 
tained by including the XLF data, which provide such a 
constraint, is M v < 0.33 eV (corresponding to the gold con- 
tours in Figure 1), nearly a factor of two improvement. 

We note that when the results excluding the XLF data 
are combined with the Gaussian prior as = 0.82 ± 0.05, 
based on the XLF results for a flat ACDM cosmology with 
n s — 0.95 (Paper I), the limits on M v obtained are virtually 
identical to those from the full combination of data, M v < 
0.33 eV. From this we conclude that the XLF data currently 
contribute to the bound on M v only through their ability 
to constrain as', the effect of non-zero neutrino mass on the 
(time-dependent) growth of structure is not detected in the 
present data. 

5.2 Extended models 

Next, we consider how the neutrino mass limits are affected 
when additional degrees of freedom are introduced in the 
cosmological model. Several possibilities are worth consid- 
ering here. To begin with, the model for dark energy may 
be generalized; the simplest ways to accomplish this by in- 
troducing a single additional free parameter are through the 
spatial curvature (i.e. by not linking the dark energy den- 
sity directly to the matter density) and through the equation 
of state, retaining the assumption of spatial flatness in the 
latter case. The flat ACDM case corresponds to effective cur- 
vature density = and equation of state w = — 1. Since 
current data do not support departures from the flat ACDM 
model either through flk / or w / -1, we consider these 
cases separately rather than introducing both parameters si- 
multaneously. Secondly, we consider the effect of primordial 
tensor modes, parametrized through the tensor-to-scalar ra- 
tio, r. We additionally marginalize over the tensor spectral 



index, n t , rather than assuming a particular model of infla- 
tion where n% would be linked to r via a consistency relation. 
Finally, we consider the effect of varying N c g on the bounds 
obtained for M v , leaving a discussion of the constraints on 
N e g itself for Section 6. 

The effects of adding each of these degrees of freedom 
individually to the ACDM+M„ model are displayed in Fig- 
ure 2 (95.4 per cent confidence regions only). The left panel 
shows results from the combination of CMB, SNIa, BAO and 
/gas data. Here it is clear that degeneracies exist between M v 
and some of the additional parameters. In particular, the 
presence of primordial tensor modes (Mv < 0.92 eV) or cur- 
vature (M v < 1.44 eV) significantly degrades the limits on 
neutrino mass relative to the simpler model (M v < 0.61 eV) 
by lengthening the major axis of the confidence region. In 
constrast, while N c s has a small effect on the M v limit 
(Mv < 0.67 eV), it significantly degrades the constraints 
on other cosmological parameters, including as, resulting 
in the much thicker confidence region in the figure (see also 
Section 6.1). The equation of state of dark energy is an in- 
termediate case which increases the bound on the neutrino 
mass somewhat (M v < 0.75 eV), while also degrading the 
constraints on other parameters slightly. Marginalized con- 
fidence intervals on parameters of interest for these cases 
can be found in Table 1. 

The right panel of Figure 2 shows the equivalent tests 
when the XLF data are included in the analysis. In this 
case, the weakest constraint occurs when N e a is free (M v < 
0.54 eV); this is sensible, since it had the smallest correla- 
tion with as a priori (left panel). In the other cases shown, 
limits on M v are improved by a factor of 2-3. In the most 
general model considered, where curvature, tensors and ad- 
ditional relativistic species are all marginalized over, the full 
combination of data yields M v < 1.31 eV, compared with 
M v < 2.34 eV when the XLF is not used (Table 1). 
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Sm„ (eV) Sm„ (eV) 

Figure 2. [preprint note: a monochrome version of this figure appears after the references, as Figure 6] Joint 95.4 per cent confidence 
regions for M y and <rg for various cosmological models. Yellow contours correspond to the basic ACDM+M„ model, blue contours are 
marginalized over Q.^, red contours over w, green over r and nt, and purple over N e g. The left panel shows constraints obtained from 
the combination of CMB, /gas, SNIa and BAO data; the right panel shows results that include the XLF in addition to those data. No 
external prior on Ho is used. 



Table 1. Constraints on the species-summed neutrino mass and other model parameters, including conservative systematic allowances. 
Note that limits on M v and r are listed at 95.4 per cent confidence, while all others are 68.3 per cent confidence. Parameters with a 
single value listed were fixed at that value. Constraints on Q m do not appear in the table, as they were extremely similar for all models, 
consistent with f2 m = 0.26 ± 0.015. The first column indicates whether the XLF data were included in the fit; the combination of CMB, 
/gas, SNIa and BAO data is used in all cases. The results in this table do not use an external prior on Hq. The use of such a prior, in 
addition to the XLF data, greatly improves constraints on both N e s and M„ when N e g is a free parameter (see Table 2). 
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Figure 3. Joint 68.3 and 95.4 per cent confidence regions on 
Q m h 2 and N e fi from the combination of CMB, / ga s, SNIa and 
BAO data with a direct measurement of Ho (blue), and the same 
data with the addition of the XLF (gold). Also shown are the 
constraints from the former combination of data, but without the 
prior on Ho (green); this illustrates the strong sensitivity of the 
N e ff results to the Hubble parameter. The dotted, horizontal line 
indicates the standard value of iV e ff , 3.046. 
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Figure 4. Joint constraints on Q^/Um and Ho when jV e g is free. 
Green contours show results obtained from the combination of 
CMB, /g as , SNIa and BAO data; despite the inclusion of SNIa 
and BAO data, the strong degeneracy in this plane that occurs in 
both / gas and CMB analyses (when N e ff is free) remains evident. 
The addition of a measurement of the Hubble parameter at the 5 
per cent level (blue contours) significantly improves the results. 



6 CONSTRAINTS ON N cB 
6.1 Simple models 

Joint constraints on N e g and Sl m h 2 are shown in Figure 3 
for a flat ACDM cosmology. The strong correlation in this 
plane is due to the fact that the CMB data constrain di- 



rectly 



which is a degenerate combination of N B ft and 



Q m h 2 , as discussed previously in Section 2. The combina- 
tion of data that we use places tight constraints on fi m , so 
significant improvement can be obtained by incorporating a 
direct measurement of the Hubble constant; throughout this 
section we use a Gaussian prior, h = 0.742 ± 0.036, based 
on the results of Riess et al. (2009). The constraints thus 
obtained are listed in Table 2. 

We note that, ordinarily, the combination of / gas and 
CMB data provides a tight constraint on Ho (Allen et al. 
2008); however, this is not the case when N B ff is free. The 
reason why can be seen in Figure 4. The green contours 
in the figure show constraints obtained from combining the 
CMB, / g as, SNIa and BAO data (without a prior on H ). A 
correlation in this plane occurs naturally in the / gas analy- 
sis, since X-ray observations of clusters measure a degenerate 
combination of the cosmic baryon fraction, /fi m , and dis- 
tance. When N c fc is fixed, CMB data place tight constraints 
on the baryon fraction, though not on the Hubble parame- 
ter, and so the combination of CMB and / gas data produces 
tight constraints on both Ho and fib/^m- When N e s is free, 
however, the CMB constraints in this plane are degenerate 
along nearly the same axis as the / gas constraints, 11 and so 
the inclusion of additional, independent distance measure- 
ments is necessary to place a constraint on Ho- The com- 



11 Respectively, the quantities f2 m /i 2 and f2 m /i 15 are degenerate 
in the analysis of CMB and / gas data. 
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Figure 5. Joint 68.3 and 95.4 per cent confidence regions on 
Q m h 2 and erg from the combination of CMB, / gas , SNIa and BAO 
data with a direct measurement of Ho (blue), and the same data 
with the addition of the XLF (gold) when JV e g is free. The tight 
constraint on as provided by the XLF data improves the deter- 
mination of f2 m /i 2 , which translates to an improved constraint on 
N cS (Figure 3). 



bination with a direct measurement of Ho (blue contours) 
significantly improves matters. 

As Figure 3 shows, the addition of the XLF data im- 
proves the constraint on N c s somewhat, from N c s = 3.61q \ 
to iV e ff = 3.4^Q g (68.3 per cent confidence). The mechanism 
for this improvement is a degeneracy between Qmh and as, 
shown in Figure 5, which the XLF data reduce through their 
constraint on as- 
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Table 2. Constraints on the neutrino effective number and other model parameters. Note that limits on M v and r are listed at 95.4 per 
cent confidence, while all others are 68.3 per cent confidence. Parameters with a single value listed were fixed at that value. The first two 
columns respectively indicate whether the prior on Ho and the XLF data were included in the fit; the combination of CMB, / ga s, SNIa 
and BAO data is used in all cases. 
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6.2 Extended models 

The addition of nuisance parameters in the form of curva- 
ture, tensors, and non-zero neutrino mass tends to weaken 
the constraints on N e g and Ho along their primary degen- 
eracy axis, similarly to what was observed in Section 5.2. 
Table 2 lists the constraints obtained with and without the 
Ho prior when various nuisance parameters are marginal- 
ized over, as well as when the XLF data are included in the 
fit. As the table shows, the inclusion of the prior on Ho ef- 
fectively eliminates the degeneracies between JV e fr and the 
nuisance parameters; in particular, the constraint on jV c ff 
obtained when marginalizing over curvature, tensors and 
neutrino mass simultaneously, N c g = 3.7+ u 'g, is very sim- 
ilar to those obtained when the nuisance parameters are 
fixed, N e ff = 3.6l"Qg. With so much extra freedom in the 
model, the inclusion of the XLF data produces only a small 
improvement in N e g in that case. 



7 PROSPECTS FOR IMPROVEMENT 

An interesting question is where future improvements in the 
determination of M v and N e g from cosmological data will 
come from. One simple way of addressing this is to see what 
other cosmological parameters are most correlated with the 
parameters of interest in the current results. We consider the 
simple cases of a ACDM model where either M„ or Neg is 



additionally free, as well as the more general case with both 
parameters free in addition to curvature and tensor modes. 

7.1 Improvements on M v 

When only M„ is additionally free, it remains most degen- 
erate with cr$, despite the tight constraints placed on this 
parameter by current data, with a correlation coefficient of 
p — —0.63. In the more general model (with M v , N e g, Q,^, r 
and nt free) , the degeneracy with as is also relatively impor- 
tant (p = —0.51), but degeneracies with fl m h 2 (p — 0.60) 
and Slk (p — 0.48) appear at approximately the same level. 

The situation will change, however, with the avail- 
ability of improved CMB data from Planck, and ul- 
timately from future high-resolution CMB polarization 
and lensing observations. These new data will tighten 
constraints on a number of parameters that contribute 
to the width of the degenerate region in the M v -os 
plane, 12 (Figures 1 and 2; Kaplinghat, Knox, & Song 2003: 
Bashinsky & Seljak 2004; Colombo, Pierpaoli, & Pritchard 
2009; de Putter, Zahn, & Linder 2009) resulting in a 
stronger correlation of M v with erg. 

12 In addition to the curvature and tensor parameters discussed 
in this paper, standard cosmological parameters such as Qmh 2 , T, 
n B and the high-rcdshift power spectrum amplitude all contribute 
to the width. 
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Using a simulated Planck data set produced using the 
FuturCMB code of Perotto et al. (2006), we project that 
such data, combined with a 2 per cent constraint on as, 
will produce a limit M v < 0.18 eV for minimal neutrino 
mass in the normal hierarchy (M v = 0.056 eV). Alterna- 
tively, M u = would be ruled out at approximately la 
significance for minimal neutrino mass in the inverted hi- 
erarchy (M v — 0.095 eV). A constraint on as at this level 
may be possible in the near term from galaxy cluster ob- 
servations, primarily by better understanding cluster mass 
measurements (e.g. Mahdavi et al. 2008; Zhang et al. 2010; 
von der Linden et al., in preparation), as well as from 
improved measurements of the CMB power spectrum at 
multipoles ~ 2000 (Lueker et al. 2009; Fowler et al. 2010; 
Komatsu et al. 2010). 

As growth of structure data become available at higher 
redshifts, the time-dependent effects of neutrino mass on 
structure formation will also provide a powerful and more 
direct constraint on M v . Future cluster surveys may pro- 
vide limits an order of magnitude better than current results 
(Wang et al. 2005). Eventually, other cosmological observ- 
ables, e.g. the cross-correlation of Lyman-alpha absorption 
with CMB convergence (Vallinotto et al. 2009), may extend 
measurements of structure formation to even higher red- 
shifts {z ~ 2). 



7.2 Improvements on N e s 

Unsurprisingly, the correlation of N c g with fi m ft 2 is dom- 
inant in both the simple (ACDM+A^ff) and more general 
(M v , JVcff, fik, r and n t free) models, respectively with 
p = 0.94 and 0.87, although we note that in the simple case 
there is also a strong degeneracy with as (p = 0.73). Thus, 
improved measurements of the Hubble parameter and/or 
mean matter density will be very significant for cosmologi- 
cal constraints on JV e ff ■ 

In particular, we consider the possibility of a direct 
constraint on Ho at the 2 per cent level. Such a con- 
straint could come from several avenues, including Spitzer 
(and ultimately the James Webb Space Telescope) obser- 
vations of Cepheids and SNIa (Freedman 2009; Riess et al. 
2009), gravitational lensing time delays (e.g. Saha et al. 
2006; Oguri 2007; Coe & Moustakas 2009; Dobke et al. 
2009), or the comparison of the Compton-j/ parameter in- 
ferred from Sunyaev-Zel'dovich and X-ray observations of 
clusters (Schmidt, Allen, & Fabian 2004; Bonamente et al. 
2006; Rapetti, Allen, & Mantz 2008). The latter two nat- 
urally complement efforts to extend measurements of the 
growth of structure using clusters to larger samples and 
higher redshifts. Even with no other improvements in cosmo- 
logical data, such a measurement (we take h = 0.72±0.0144, 
for concreteness) would provide a constraint N c g = 3.04 ± 
0.35 in the simple ACDM+A cff model, or N cB = 3.14 ±0.48 
in the general model (parameters listed above; 68.3 per cent 
confidence) . Since N e fi and M v are degenerate in this model, 
the upper limit on M„ would also be improved, from 0.70 to 
0.55 eV. 



8 CONCLUSION 

We have applied measurements of the growth of massive 
galaxy clusters (detailed in Papers I and II), in combination 
with other cosmological data, to the problem of constrain- 
ing the species-summed neutrino mass, M v , and effective 
number, N e s- 

Our results show that a robust measurement of erg from 
clusters significantly improves limits on M v from current 
data, and reduces their sensitivity to assumptions about the 
cosmological model. In a simple ACDM+M„ cosmology, the 
addition of the XLF data improves the 95.4 per cent confi- 
dence upper limit on M v to 0.33 eV, compared with 0.61 eV 
from the combination of CMB, / ga s, SNIa and BAO data. 
In a more general model, marginalized over spatial curva- 
ture, primordial tensors and the effective number of neutri- 
nos, and incorporating a prior on the Hubble parameter, the 
XLF data improve the limit from 1.75 eV to 0.7 eV. The re- 
sults indicate that this improvement is due entirely to the 
ability of the cluster data to constrain as (at 2 = 0); while 
measurements of the (z-dependent) growth of structure in 
principle contain additional information about M v , current 
data are not sufficient to exploit it. 

The cluster data additionally improve constraints on the 
effective number of relativistic species, from N e ff = 3.6lg \ 
to N eS = 3Ato 6 s in a simple ACDM+A^a model (68.3 per 
cent confidence). For the more general model, including cur- 
vature, tensors and neutrino mass, only a marginal improve- 
ment is observed. 

The results obtained here are compatible with other re- 
cent estimates based on galaxy clusters (Reid et al. 2010; 
Vikhlinin et al. 2009), reflecting the good agreement in re- 
cent as constraints based on X-ray and optically selected 
clusters (Paper I; Henry et al. 2009; Rozo et al. 2010). Al- 
though current cosmological data are not sufficient to rule 
out the existence of sterile neutrino species or distinguish 
between the normal and inverted mass hierarchies, they con- 
tinue to provide some of the tightest constraints available, 
in particular on the mass scale. 

We also consider the prospects for further improvement, 
finding that a few per cent level measurement of the Hubble 
parameter will significantly improve the constraints on mod- 
els where N e g is free. A similarly improved determination of 
as, combined with improved CMB data from Planck, will 
further tighten limits on the neutrino mass scale, 13 perhaps 
providing the first hints of non-zero mass from cosmological 
data. More precise measurements of the growth of structure 
at late times, and the extension of such data to higher red- 
shifts, will provide a powerful, new mechanism to constrain 
neutrino mass. 



13 We note that the 7-year WMAP results, which were released 
while this work was in revision, already offer some improvement. 
Specifically, Komatsu et al. (2010) find M„ < 0.58 eV by combin- 
ing 7-year WMAP data with the Riess et al. (2009) Ho prior and 
the BAO results of Percival et al. (2010), a noticeable improve- 
ment over the limit M v < 0.66 eV obtained by Sekiguchi et al. 
(2010) from the same auxiliary data sets combined with WMAP5. 
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Figure 6. (Black and white version of Figure 2.) Joint 95.4 per cent confidence regions for M u and erg for various cosmological models. 
Solid contours correspond to the basic ACDM+Mj, model, dotted contours are marginalized over f2 k , dashed contours over w, dot-dashed 
over r and rit, and 3-dot-dashed over iV e jf. The left panel shows constraints obtained from the combination of CMB, /gas, SNIa and 
BAO data; the right panel shows results that include the XLF in addition to those data. No external prior on Ho is used. 
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